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The ground state energy of a collection of n pions and m kaons with short range interactions 
is calculated for a volume with finite spatial extent L and periodic boundary conditions. This 
calculation is accomplished to order in the large volume expansion. With this result one can 
extract the various two- and three-body interactions between pions and kaons from lattice QCD 
data. 
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o ■ 

■ Recently, efforts have been made using lattice QCD to extract both the Tr+Tr"*" and the K^K^ scattering lengths as 
well as the three-body interactions 7r+7r+7r+ and X+if +i4r+[l|, [3, 4|. 



> 

^ ' lattice QCD to rigorously determine the nuclear equation of state (NEOS) directly from the underlying theory of 



J|. These efforts have made steps toward using 



QCD. A knowledge of the NEOS would be of great importance to many research areas, including helping to predict 



7^ 

' the evolution of supernovae. In addition, an understanding of multi-pion and multi-kaon systems provide insight into 



strongly interacting boson gases. References l3, l3|, have made use of a method involving the volume dependence 
~r ' of the energy spectrum (below inelastic thresholds) of two hadrons as a function of their scattering lengthP, [g], where 

Oh: 

D _ lattice QCD is used to calculate the energy of multi-hadron states and this result is then applied to the calculated 



dependence of the finite volume ground state energy shift to determine interaction strengths. The ground state ener, 

of n identical bosons with short range interactions is calculated in the large volume expansion to 0{L^'^) in Ref. [7|, 

' providing a way to extract the three-bodyinteraction strength between identical bosons which enter at 0{L^^), and 

[ building upon the earlier work of Refs. P, 0, [l^ [ll|, [l^l- In this work we extend the calculation to determine the 

^-H ' ground state energy shift in finite volume of n pions and m kaons to 0(L^^) in the large volume expansion using a 
^' Pfl 

(yQ multispecies extension of the techniques used in Refs. [3, |8[. This result will allow for the systematic extraction of 

o . 

• • • ttK scattering lengths and the three-body interactions ttttK and ttKK, as well as provide a different way to extract 



TTTT and KK scattering lengths by using mixed pion-kaon systems in lattice calculations. 

iival( 
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■ To calculate the ground state energy there are multiple (and entirely equivalent) methods that one can use. One 



such method uses a resolvent function of the system Hamiltonian defined by[ 

Fiz) = (0|^|0) 



z-Eo- r{z) 
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with 

riz)^{0\Vp^(^-^V^ |0) (2) 

where |0) is the ground-state of the system and Qo = 1 — |0)(0| is a projection operator discussed below. Here H is 
the full Hamiltonian of the system while Ha is the free Hamiltonian. In eqn. ([1]), when F{z) takes the form found 
in the first line, a perturbative expansion leads to a pole near z — Eq. However, projecting out that portion of the 
propagator containing the pole using the operator Qo leads to the function r{z) in the second line which is smooth in 
the neighborhood of Eq and which can be perturbatively expanded in powers of the potential V. To find the ground 
— .... „„e .... .„ ... ...... ... po.. . . T.. .e™ .„ .M. 
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jl dzi 

AE = ro{Eo) + ro{Eo)n{Eo) + ro{Eo)rl{Eo) + rl{Eo)r2{Eo) + ... . (3) 



The method of pseudo-potentials can be employed to determine r{z) above, or one can use the different but 
equivalent approach of nonrelativistic time- independent perturbation theory (NRPT). The potential used in this 
calculation is: 

E ~ rj)5^^\n - ru) + ^ V3,^7.K6^^\n - r.j)5^^\u - rk)+ 

E V3,.KKS^^\n-rj)S^^Hn-rk)+ J2 V3,KKKS^'Hn-rj)S^'\n-rk) (4) 
iev^,j<keVK i<j<keVK 

where the first line is the two-body interaction potential while the remaining lines constitute the various three-body 
potentials. The sets Vtt and Vk refer to the set of all pions and kaons in the system, respectively. Due to the 
requirement that bound states do not form, the interactions between pions and kaons studied are restricted to be 
repulsive. Hence all tt and K fields referenced are shorthand for tt^ and . Because the particles are nonrelativistic, 
interactions such as KK — > tttt or tttt — > KK are not considered. 

The x's can be related to the scattering amplitude and for small external momentum be expanded using effective 
range theory according to x = — ffp^^ tan S{p) ^if (^^ + h^f^P^ + ■ ■ ■ where a is the s-wave scattering length 
and To is the effective range. The modern language and description of this approach is pionless effective field theory 
or In this framework, the two-body interaction x corresponds to the expansion of the 2 — )■ 2 



contact interaction in the Lagrangianjl4|. The language of EFT('^) possesses a more natural way of dealing with 
issues of renormalization and power counting as the choice of a particular subtraction scheme leads to order-by-order 
renormalization of loop divergences. 

The contact potential from eqn. (|4]) can be reexpressed in terms of a Hamiltonian in NRPT, as described in ref. 
Using this method, each term in the large volume expansion of the function r{z) can be expressed as a sum of 



diagrams in the perturbation theory. Vertices in these diagrams are given by the x 77 shown above, while the 
propagators are given by the expectation of the free theory resolvent y^) over intermediate states. The presence 
of the projection operator Qo removes from consideration any diagrams where the intermediate state expectation 
gives (Tio) = 0. Loop divergences are regulated with dimensional regularization in our calculation. The divergences 
arising from two-loop diagrams will necessitate the inclusion of three-body interaction strengths as a function of the 
renormalization scale ^. In finite volume the integrals normally associated with the continuum expression for each 
diagram will become sums because the momenta will be restricted to the possible values p = 27rn/i for Uj e Z. 

It is simple to show that generating all of the diagrams and summing their contributions leads to the following 
result for the ground state energy shift: 



A£'o('T., TO, L) — Eo{n,m, L) — n rriT, — m rriK 

= ET^{n, L) + EK{m, L) + ETrK{n,m, L) 



(5) 



with {iriTrK — mT^mK / (niTt + mK)) 



E,{n,L) = 



n 

3 



(X^ + (2n - 7)XJ^ + (5n2 - 41n + 63)/C) 



(6) 



E^K{n,m,L) = 



1 - 



\ttL ) 



I 



a-irK 

ttL 



J 



, , / 1 2 
— (n- 1) + — 



O-K / N / 1 2 
+ — l^(m- 1) + 




nm(n - l)773_^^if (L) nTn{m - 1)i]3^^kk{L) 



2X6 



2L6 



with the scattering lengths given in terms of the parameters Oi = and the effective ranges by 



(7) 





- a. ^3 










7=3/2 


- a^K ^3 


(8) 



The four volume dependent (but renormalization scale independent) quantities from the three-body interactions are 



defined by (y = mT^/mK)'- 



mAL) = ^3,,(^) + ^l!I^(3V3-47r)log(/.L)-^(2g[l,l]+r[l,l]) 



4a4 2 



4=0 p=7r,ft A/gA/2 



and 



M = {Q{l,y),Q{y,l),R{y,l),R{l/y,l/y)} 

N2 - {g(l,l/y),g(y,y),i?(y,y),i?(l,l/2/)}. (10) 

The functions Q, R, q, and r are defined in the appendix along with the coefficients The finite parts of Q(a, b) 
and R{a,b) are scheme dependent quantities where changes in the value will be absorbed by ri3{fi). The numerical 
values for the MS scheme are given in the appendix. However, the fj^, are not scheme dependent and this is the 
quantity that will be determined during a lattice calculation. Furthermore, note that the three-body interactions in 
the ttttK and t:KK cases depend on the ttA' mass ratio. One can immediately see that this is necessary if one takes 
the limiting case where either the pion or the kaon become infinitely heavy, where the heavy particle decouples from 
the theory and hence all cross species interactions must go to zero. Finally, in the limit of n ^ or to — ^ this result 
simplifies to the previously determined n-boson case while the limit to — )■ n with mx to^t and all interactions set 
to be equal it simplifies to the 2n-boson case. The result given agrees with previously calculated results in the single 
species limit 

In this work we have calculated the ground state energy shift of the mixed species system of n pions and to kaons. 
Using this result a rigorous connection can be made between Euclidean-space lattice QCD calculations of mixed pion- 
kaon systems and Minkowski-space multibody interaction strengths. In Ref. 3] the isospin and strangeness chemical 
potentials are analyzed with lattice QCD at values relevant to dense nuclear matter for purely kaon condensates. 
With the results presented in this paper such an analysis can be done for mixed pion-kaon systems, allowing lattice 
calculations of the t:K scattering length and three-body interaction strengths ttttK and ■nKK. Such extractions are 
necessary for examinations of the nuclear equation of state for values of the isospin and strangeness chemical potentials 
where it is not energetically favorable to form either pure pion or pure kaon condensates, but where possibly mixed 
systems of pions and kaons exist. 
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APPENDIX 



The coefficients fi are given by (where y = rriT^/mK and rriT^K = m^mK/iin^ + tok)): 



Jo 



fi 



ft 



1 



(1 + m + 



2fl + ^) (l_n) 

- rriK , 



= mT,K- 



mTrfriK 



-(1-n) 



IJ,K rXJ,TT 



(A-1) 

(A-2) 

(A-3) 
(A-4) 



Jo 



fl 



cK,K 



ml + m 



1 / 

- 4n + 4m — 3mn — 6 + 

2 V m^rriK 



— {m + n — mn — 1) 



= -8(n-l)(m-l) 1 + 



m: 



mT^rriK 



(n-1) 



/2 — 



(1 - n) [m|-(13n- 45) - ml +TO^m|^(14n - 39) + 5mlmKin - 3)] 



m, TUj^ ^ rriK) 



(n- l)(n- 2) 



(A-5) 
(A-6) 
(A-7) 

(A-8) 
(A-9) 



j-a(l,J/),7r _ j.K(i/,l),7r _ 

Jo ~ /o 



FC(l,y),7r 



Q(i,i/i/),if 



e(i,3/),-ff 



7?.(2/,l),7r 



m-rr 



a(i,i/3/),T 



TC(l,l/!/),7r 



(A-10) 
(A-11) 
(A-12) 
(A-13) 
(A-14) 



Jo 
J2 

fQ{y,v),K 

J i 



fQ{yA),K _ fQ{y,l),K 
JO — Jl 



f. 



a{yA),K 



yQ(l/,S/),'r _ Q 



(A-15) 
(A-16) 
(A-17) 
(A-18) 



f 

Jo 

My,y),K 
Jl 

n(i/yMy),'K ^ jn{y,y),K 



rit{y,y),7T _ .TC(2/,2/),7r _ rn(y,y),K _ ,n(y,y),K _ ,n(l/y,l/y),K 
Jl — J2 ~ Jo ~ Jl — Ji 



= 



m 



K 



(A-19) 
(A-20) 
fA-21) 



The integer sums that contribute to 0{L ^) inchide: 



K. 
Q{a,b) 

n{a,b) 



1 \ 

I = Urn V - 47rA, = -8.9136329 

J = ^ 16.532316 

y- 



= 8.4019240 



V V L 



E 



^ Itt-I 



3^ |n|4 |nP + a|m|2 + 6(n + a + 6 



(A-22) 
(A-23) 
(A-24) 
(A-25) 

(A-26) 



where the Q and TZ sums need to be regulated (defined) using dimensional regulation and in the TZ sum the subtracted 
integral removes a nested divergence that would otherwise affect the one-loop scattering amplitude. One can calculate 



1 



Q(a,6) = Qi:{a,b) + QDRia,b)-QA{a,b) 

A. , X 27r3 / b 

= Q(a,b,ii)-\ — sm — , , 

b \^{a + b){l + b)J e 

n{a,b) = n^{a,h)+nDR{a,b)-nA{a,b) 

, , 27r^\/a + b + abl 
= iZ[a, b, /i) — 



(A-27) 



(A-28) 



(a + 6)2 e 

where the cutoff dependencies of the sums and T?.^ are canceled by the cutoff dependence of the Qa and TZa 
integrals. The dimensionally regularized integrals in turn counter the remaining portions of the cutoff integrals, 
leaving one with the functions Q and TZ found in eqn. These are given by 

b 



Qia,b,n) = q{a,b) + —sin ^ == 
b \^{a + b){l + b)^ 



logifxL) 



7^(a,6, /i) = r(a,b) 



8n^y/a + b + ab 



log(^L) 



(A-29) 
(A-30) 



(a + 5)2 

where in Tables HI and HIl the numerical portions of q{a, b) and r{a,b) are explicitly evaluated for values of the arguments 
a and b relevant to the calculation. 

TABLE I: Numerical evaluation of the function q(a,b). 



y = — 


g(i,y) 


9(1, 1/y) 


g(y,i) 


<i{y,y) 


0.2 


-167.96763 


-34.22227 


-94.42134 


-113.71723 


0.2827 


-157.34173 


-45.55930 


-95.75651 


-127.15818 


0.4 


-144.45515 


-59.06850 


-97.36636 


-129.53719 


0.6 


-126.85940 


-77.19563 


-99.49508 


-121.84065 


0.8 


-113.15457 


-91.10749 


-101.03837 


-111.68392 


1 


-102.15598 


-102.15598 


-102.15598 


-102.15598 
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TABLE II: Numerical evaluation of the function r(a,b). 



y= UhL. 






'"(1/2/, 1/2/) 


r{y,v) 


0.2 


49.695247 


-0.358353 


-0.141373 


517.146558 


0.2827 


44.958613 


0.993223 


0.391640 


265.289351 


0.4 


38.951624 


3.803394 


1.738871 


133.442877 


0.6 


30.564875 


9.279473 


5.616888 


58.127093 


0.8 


24.132625 


14.530094 


11.449288 


31.463511 


1 


10.186903 


10.186903 


10.186903 


19.186003 
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